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Abstract 

Selectron decays constitute a source of 100% polarised neutralinos, whose 
helicity is fixed by the charge and "chirality" of the decaying selectron. In SUSY 
^ ■ scenarios where the second neutralino X2 nas three-body decays, the cascade 

decay e~i — * &X2 — ► provides a clean place to study CP violation in the 

neutralino sector, through the analysis of CP-violating asymmetries involving 
qq ' the X2 s P m * an( i the momenta of the two muons. We show that a CP-violating 

asymmetry in the triple product s ■ x could be observable at a 800 GeV 
CN ' linear collider provided the gaugino mass Mi has a large phase at the electroweak 

^ ' scale. 
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Supersymmetric theories [1-3] are perhaps the best motivated extensions of the stan- 
dard model (SM). If supersymmetry (SUSY) is realised in nature it must be broken, 
possibly at a relatively low energy scale. SUSY-breaking interactions are usually as- 
sumed to be flavour-blind, at least approximately. Otherwise, they would lead to unac- 
ceptable rates for flavour-changing neutral currents. Analogously, large phases in the 
SUSY-breaking terms of the Lagrangian give supersymmetric contributions to electric 
dipole moments (EDMs) far above present limits. For a relatively light SUSY spectrum, 
the only solutions to overcome this problem are, either assume that all SUSY-breaking 
parameters have very small phases (or are real) , or to arrange cancellations among the 
different contributions to EDMs to satisfy experimental limits [4-6]. 

Although the assumption that SUSY-breaking terms are real is conceptually simpler 
and more attractive, the possibility of large phases and apparently "fine-tuned" internal 
cancellations must not be discarded. Indeed, in quantum field theory the Lagrangian 
parameters are complex in general, unless there is some argument (e.g. hermiticity of 
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the Lagrangian or some symmetry) requiring them to be real. Thus, setting the phases 
of SUSY-breaking terms to zero "by hand" may also be regarded as fine tuning, in the 
absence of a symmetry principle to explain why these terms must be (approximately) 
real. The possibility of large SUSY CP-violating phases makes compulsory to explore 
their effects in phenomenology, in order to find out their presence and determine their 
magnitude. 

In this Letter we are interested in the direct observation of CP violation due to 
supersymmetric CP phases, rather than in their effect in CP-conserving quantities 
such as cross sections and decay widths (see for instance Refs. [7-10]). We restrict 
ourselves to the minimal supersymmetric standard model and focus on the neutralino 
sector, studying CP asymmetries in the cascade decay ef — > e ± x\ — > e ± xV + / i • The 
X2 produced in decays are 100% polarised, having negative helicity in — > e~X2 
and positive helicity in e\ — > e + x°- Having perfect x\ polarisation is obviously a 
great advantage for the study of angular correlations involving the x\ s P m [H], an d in 
particular for the study of CP asymmetries. Selectrons will be discovered at LHC [12] 
if they have masses of a few hundred GeV, but the large backgrounds present make it 
impossible a detailed study of their properties, which must be carried out at a e + e~ 
collider. Therefore, as source for left selectrons we consider e^e^ and e^ei, production 
in the processes 

e + e" -> e\el -> e + x\e~xl ~* e 4 X? e ~ xV + >~ , 

e + e~ -> e+e~ -> e + x\e~xl -»• e + X?e~xW~ , 

e + e~ -> e R e\ -> e~x?e + x° -»• e^X?e + x?/i + /i" . (1) 

In these processes all final state momenta can be determined, and the selectron and 
neutralino rest frames can be reconstructed [13], allowing the study of CP asymmetries 
involving the x\ spin s and the momenta of ji~ in the X2 res t frame. We do not 
consider e^e^ production with one selectron decaying to ex\i because the branching 
ratio of this decay is very small in general. We discuss SUSY scenarios where the second 
neutralino has three-body decays, in which case it is possible to have a CP asymmetry 
in the triple product s • (p M - x p^+) of order 0.1. In scenarios with two-body decays 
X2 ~^ x\Z the asymmetry is of order 0.02, and when two-body decays to sfermions 
dominate it is negligible. We consider a CM energy Eqm = 800 GeV, as proposed for 
a TESLA upgrade. At this CM energy, the cross sections for production are higher 
than at 500 GeV due to the smaller destructive interference between s- and t-channel 
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contributions and the larger phase space available. 



2 Decay of X2 an d CP asymmetries 

The decays of the second neutralino X2 ~^ Xiff are mediated by the Feynman diagrams 
in Fig. [TJ (The diagrams in Fig. [TJb with neutral scalars are negligible except for 
/ = t,b,r because they are proportional to the fermion Yukawa couplings.) In this 
work we are interested in the final state with ff = in which the energies, 

momenta and charge of both particles can be measured. The decay channel with 
// = e + e~ shares these properties, but the multiplicity of electrons in the final state 
makes it difficult to identify the ones resulting from the X2 decay. The bb final state 
is also interesting for its large branching ratio, but the b tagging efficiency reduces the 
signal to the same cross section of the channel. 




(a) (b) 




(c) (d) 

Figure 1: Feynman diagrams for the decay X2 ~^ Xiffi mediated by Z bosons (a), 
neutral Higgs bosons (b), and scalar fermions (c and d). 

We consider SUSY scenarios where the second neutralino has three-body decays. 
This happens when all sleptons are heavier than x\ an d m ^ < M z + m^o . For defi- 
niteness, we choose a scenario similar to SPSla in Ref. [15] but with a heavier sfermion 
spectrum and complex phases in Mi and \i. In this scenario Xi an d X2 are g a ugino- 
like. The low-energy parameters (at the scale Mz) most relevant for our analysis are 
collected in Table [TJ For <pi = (p^ = 0, these parameters approximately correspond 
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to m V2 = 250 GeV, = m L = m Hi = 200 GeV, A E = -200 GeV at the unifi- 
cation scale, and tan/3 = 10. We use SPheno [14] to calculate sparticle masses and 
mixings, as well as some decay widths. Neutralino masses slightly depend on <pi and 
M . For 0i = M = they are m^o = 99 GeV, m^o = 178 GeV, m^o = 384 GeV, 
m^o = 400 GeV, m-- = 177 GeV, and for other values of </>i they differ up to ±2 
GeV. The relevant branching ratios (taking 0i = M = 0) are Br(e^ — > ex?) = 0.18, 
Br(e fl - ex?) = 0.998, Br(g L - e#j) = 0.30, Br(x° - X> + /0 = 0.035. 



Parameter 


Value 


Mi 


102.0 e^ 1 


M 2 


192.0 


A* 


377.5 e**" 


tan/5 


10 




224.0 




264.5 




252.4 



Table 1: Low-energy parameters (at the scale M z ) for the SUSY scenario used. The 
dimensionful parameters are in GeV. 

CP violation in the neutralino mixing matrix arises when the phases of Mi and/or 
/j are different from 0, 7r. These phases generally lead to large supersymmetrc con- 
tributions to EDMs. If the squark spectrum (which does not play any role in our 
analysis) is heavy enough, experimental limits on the neutron and Mercury EDMs are 
satisfied. On the other hand, for the values for selectron masses under consideration, 
the experimental bound on the electron EDM d e imposes a severe constraint on the 
allowed region in the (0i,0 M ) plane. Using the expressions for d e in Ref. [16], we find 
that for each 0i between and 2n there exist two narrow intervals for M (one with 
values M ~ and the other with values </>^ ~ 7r) in which the neutralino and chargino 
contributions to d e cancel, resulting in a value compatible with the experimental limit 
rf e ex P = (0.079±0.074) x 10~ 26 e cm [17]. (For instance, for fa = tt/2 we find M ~ -0.12 
or M ~ 3.21.) Therefore, in principle it is possible to have any phase </>i, though with 
a strong correlation with M . Of course, if <\>\ and M are experimentally found to be 
non-vanishing, a satisfactory explanation will be necessary for this correlation, which 
apparently would be a "fine tuning" of their values [18]. 

Let us discuss which asymmetries may be defined in the \\ decay. In the x° res t 
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frame, the decay looks as depicted in Fig. El with s the X2 s P m an d the 3-momenta in 
obvious notation. Under CP, the spin and momenta transform as: 



-P^ 



Pxl 



-Px? 



Xi ■ 



(2) 




Figure 2: Schematic picture of the xi decay in its rest frame 



Using these vectors, we can build the CP-odd quantities 

Qi = s- (p„- x P»+) > 

Q2 = a - (p,,- + ) . (3) 

Several other CP-odd quantites can also be constructed involving higher powers of the 
momenta or using the energies, e.g. we can define Q3 = — which is also odd 

under CP. For Q1-3 we construct the asymmetries 

= N(Q t > 0) - N(Q t < 0) 
1 N(Qi > 0) + N(Q t < 0) ' U 

where N denotes the number of events. Since x\ and Xi are Majorana particles, these 
asymmetries must vanish if CP is conserved. Hence, they are genuine signals of CP 
violation. We note that Q2,3 are even under "naive" time reversal T. This implies 
that v4 2)3 need the presence of absorptive CP-conserving phases in the amplitude in 
order to be nonvanishing [19]. CP-conserving phases appear in loop diagrams with on- 
shell propagators, through final state interactions or from the widths of intermediate 
unstable particles (the phases originated by particle widths are usually tiny). In x° 
decay, T-even asymmetries like ^3 can result from the interference of a dominant 
tree-level and a subleading loop diagram mediating the decay and are expected to be 
small. On the other hand Qi is T-odd, thus a relatively large asymmetry Ai can be 
generated at the tree level, without the need of an absorptive phase. Our analysis is 
centred in the asymmetry A\. 
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3 Generation of the signals 



We calculate the matrix element for the resonant processes in Eqs. (|T} using HELAS [20], 
so as to include all spin correlations and finite width effects. The relevant terms of the 
Lagrangian and conventions used can be found in Refs. [13,33]. We assume a CM 
energy of 800 GeV, with e~ polarisation P e - = —0.8 and e + polarisation P e + = 0.6, 
and an integrated luminosity of 534 fb _1 per year [21]. Beam polarisation does not have 
any effect on the asymmetries, which are defined for decays independently of the 
production mechanism, but increases the total signal cross sections. In our calculation 
we take into account the effects of initial state radiation (ISR) [22] and beamstrahlung 
[23,24]. For the latter we use the design parameters T = 0.09, N = 1.51 [21]. 1 We also 
include a beam energy spread of 1%. 

In order to simulate the calorimeter and tracking resolution, we perform a Gaussian 
smearing of the energies of electrons (e) and muons (fi) using the specifications in the 
TESLA Technical Design Report [25] 

AE e 10% „ AE* „ „ 

©1%, -=— = 0.02%^, (5) 



E e 1 

where the energies are in GeV and the two terms are added in quadrature. We apply 
"detector" cuts on transverse momenta, pr > 10 GeV, and pseudorapidities \r]\ < 2.5, 
the latter corresponding to polar angles 10° < 9 < 170°. We also reject events in which 
the leptons are not isolated, requiring a "lego-plot" separation AR = a/ Ar] 2 + A0 2 > 
0.4. We do not require specific trigger conditions, and we assume that the presence 
of charged leptons with high transverse momentum will suffice. For the Monte Carlo 
integration in 8-body phase space we use RAMB0 [26]. 

The reconstruction of the final state momenta is done with the procedure described 
in detail in Ref. [13], with some modifications. In general, it is necessary to have as 
many kinematical relations as unknown variables in order to determine the momenta 
of the undetected particles. In our case, there are 8 unknowns (the 4 components 
of the two Xi momenta) and 8 constraints. These are derived from energy and mo- 
mentum conservation (4 constraints), from the fact that the two Xi are 011 shell (two 
constraints), from the kinematics of the decay of the x° ( one constraint). The last con- 
straint comes from the hypothesis that in e + e~ collisions two particles are produced, 
either with the same mass or having a squared mass difference equal to m\ L —tu\ . The 
reconstruction is attempted for both cases, considering different effective CM energies 
1 The actual expressions for ISR and beamstrahlung used in our calculation can be found in Ref. [13]. 
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E e Q < E CM (to partially take into account ISR and beamstrahlung effects), selecting 
the one which gives reconstructed selectron masses closest to their actual values which 
can be measured in other processes [27-30]. If the event does not reasonably fit into 
any of the two possibilities, it is discarded. The identification of tifih versus tufiL 
production is successful in most cases, with a good "tagging" 80% of the time for e^L 
and 70% for tRtL- We note that these processes are topologically very similar, being 
the only difference the selectron energies in CM frame (and thus the magnitude of their 
3-momentum). In the former process, both selectron energies are 400 GeV, whereas in 
the latter they are 412 GeV for e L , 388 GeV for e R . 

The reconstruction procedure determines the momenta of the two unobserved x? ? 
identifying for e£e^ whether the selectron pair has decayed in the channel e\ — > e + Xi, 

— > e^x® ~^ e _ Xi/ i+ A t ~ or 5 on the contrary, in the channel e\ — > e + x\ — > e + X\^ + >~, 
e L — > e~Xi- The knowledge of all final state momenta, as well as the identification of 
the particles resulting from each decay, allows us to construct various mass, angular 
and energy distributions [11], and in particular the determination of the selectron and 
x\ rest frames. In each event, the x\ spin direction s can be found as follows: If the x° 
results from the decay of a , it has negative helicity, so its spin direction is s = —p, 
with p the X2 momentum in the rest frame. If the X2 results from the decay of a 
e£, it has positive helicity and its spin direction is s = p, with p the X2 momentum in 
the e\ rest frame. 

Finally, it is necessary to discuss the possible SM and SUSY backgrounds to our 
signal. The relevant SUSY backgrounds are sneutrino and x\ pair production, in the 
decay channels 

e + e" -> z>*z> e -> e + Xi e~xt -»■ e + xV~^ e_ Xi^ + ^ > 

e + e~ -> u*u e -> e + Xi v e xl -»■ e + X?e~z/ e z/ e V~ , 

e + e" -> -»■ X?e + e" XiP + P~ ■ ( 6 ) 

In sneutrino pair production the x^x^ decay channel has a larger branching ratio, 
with Br(i> — > e~xt) — 0.52, Br(%^ — > x5 e+z/ ) = 0-1) while for the the X1X2 m °de 
Br(z> -> 1/X2) = 0-22, Br(x°2 -»■ = °- 035 - We have calculated the three 

processes in the same way described for selectron pair production, and include them in 
our results. The cross sections of these backgrounds turn out to be 4 times larger than 
the signals, but can be reduced with the reconstruction of the final state momenta. The 
reconstruction method applied for the signal partially eliminates the three backgrounds. 
The first and third ones are further reduced applying a reconstruction procedure specific 
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for each case (not discussed here for brevity) and requiring that signal and background 
events do not have a kinematics similar to sneutrino or X2 P a i r production. 2 This is 
not possible for the second background in Eq. (J@J) due to the different kinematics of 
the process. A cut requiring transverse energy Ht > 200 GeV is applied as well. Other 
SUSY backgrounds like e + e~ — > XxX^ ~^ XiXi^ — > e + e~ fM + vvXiXi (involving 
several decay channels which lead to the same final state of e + e~/i + yU - plus missing 
energy and momentum) are much smaller, with cross sections smaller than 0.01 fb. 
The SM background is given by six fermion production e + e~ — > e + e~ /i + pT vv . Its 
cross section is below 0.3 fb [31] and with kinematical cuts can be eliminated more 
easily than SUSY backgrounds. 



4 Numerical results 

We first examine the possible values that this asymmetry may take in connection with 
EDM constraints. In Fig. El we show its dependence on the two phases fa, M (these 
plots are calculated with the Monte Carlo described in last section, but assuming perfect 
momenta reconstruction and particle identification, and without any kinematical nor 
detector cuts). For fa sufficiently large A\ reaches values of ±0.13, while for fa = the 
asymmetry is negligible independently of M . 3 The dependence of the cross section on 
these two phases is plotted in Fig. HJ Our approach is then as follows: for each value 
of fa, we know from the discussion in Section 01 that there is an allowed interval of M 
(which we may take with < 0.12) in which the electron EDM does not exceed the 
experimental bound. We then calculate A\ for this (fix, but taking M = 0, since the 
asymmetry is almost independent of the latter phase and the effect on the cross section 
is also rather small for j^j < 0.12. 

The cross sections for fa — of the different processes in Eqs. (I1IGI) are collected 
in Table IH before and after final state momenta reconstruction and kinematical cuts. 
All figures include the various corrections discussed in the previous section. The recon- 
struction of the signal and backgrounds allows us to reduce the latter by a factor of 20, 

2 In the X1X1 channel only the momenta of each vxi P & i r can be determined, but that is enough 
to obtain the chargino and sneutrino momenta. 

3 The main contribution to the asymmetry comes from the interference between the Z diagram in 
Fig. GJi and the diagrams with exchange in Fig.^J; andOJi. We note that the Z contribution alone 
produces an asymmetry A\ — 0.021 without the need of interference, but the asymmetry is much 
larger when the rest of diagrams are included. 
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(a) (b) 

Figure 3: Theoretical value of A\ as a function of 0i, for M = (a) and as a function 
of M , for 0i = (b). 




Figure 4: Total cross section (without ISR, beamstrahlung and beam spread effects, 
and without any kinematical cuts) as a function of (fix, for M = (a) and as a function 
of M , for 4>i = (b). 

while keeping 80% of the signal. We define the normalised quantity Q\ = s-(p M - x 
using unit vectors in the direction of the muon momenta, so that — 1 < Q\ < 1. The 
kinematical distribution of Qi for the signals is shown in Fig. El taking fa — and 
= it/2. We have normalised the e^e^ cross sections to unity, and the e^ej ones to 
1/3. In this plot we take into account ISR, beamstrahlung, beam spread and detector 
effects. The Qi distribution already shows that the asymmetry is a real dynamical 
effect and is not a fake asymmetry caused by "detector" cuts applied in phase space. 
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In our calculation we find that for 0i = 7r/2 the relative difference between the phase 
space volumes of the two hemispheres (with Qi < and Q\ > 0) is negligible, of 
3.1 x 1CT 4 for e£e^ and 6.9 x 1CT 4 for e^ej. We note that this is not the case for 
the asymmetry A 2 . Even in the CP-conserving case with <p% = a fake asymmetry 
A 2 = —0.055 is generated just by phase space cuts, and might obscure the observation 
of a real CP asymmetry, which is expected to be very small in this case. 

Before After 



elel 0.21 0.18 

e+eZ 0.061 0.046 

e+e^ 0.027 0.021 

v*Z (Xlxt) 0.99 0.039 

W (XiXi) 0.24 0.011 

Xlxl 0.11 0.0081 



Table 2: Cross sections (in fb) before and after signal reconstruction and kinematical 
cuts of the processes in Eqs. (jl|6j) . 




The difference between the distributions for the two values of the phase is smaller for 
the tRth signals. This is because the reconstruction procedure has a smaller efficiency 
for the correct identification of e^ej events, leading to a small washing-out of the 
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asymmetry. For 4>i — 7r/2, we have for e^e^ the value A\ = 0.110, while for e^e^ 
and we find A\ = 0.098, Ay = 0.099, respectively. The total asymmetry A\ as a 
function of <fii is plotted in Fig. Ek, including the backgrounds in Eqs. (JBJ), which have 
vanishing CP asymmetry, as well as ISR, beamstrahlung, beam spread and detector 
effects. The shaded region represents the statistical error for two years of running, with 
a luminosity of 534 fb _1 per year. The cross section also exhibits a strong dependence 
on 0i, as can be observed in Fig. Eh. Nevertheless, the measurement of the cross section 
is not likely to provide any information on the phase 0i, due to the multiple theoretical 
uncertainties present regarding the neutralino mixing matrix, sparticle masses, scale 
dependence, etc. 



0.45 - 




Figure 6: (a) Asymmetry A\ as a function of the phase 0i. The shaded region represents 
the statistical error for two years of running, (b) Total cross section for e + e~ — ► 
&lr&l ~^ e+e ~/ i+ A i ~XiXi as a function of 0i. In both cases the backgrounds from u*v 
and X2X2 production are included. 



5 Summary 

In this work we have shown that a CP asymmetry in selectron cascade decays — > 
e X.2 ~> e Xi/ i+ / i ~ can be observable in SUSY scenarios where X2 has three-body decays, 
if there is a large nonzero phase 0i in the gaugino mass M\. The key features for 
the observation of the CP asymmetry are: [%) the neutralinos produced from selectron 
decays are 100% polarised; [%%) all final state momenta can be accurately reconstructed 
in the processes under consideration; {%%%) this reconstruction allows us to eliminate 



11 



potentially dangerous backgrounds. 

In SUSY scenarios with three-body decays of a sizeable production can only 
take place at a TESLA upgrade with a CM energy of 800 GeV. We have selected one 
of such scenarios, with Xi an d x\ gaugino-like, and shown that the CP asymmetry in 
the triple product s • x p^+) could be up to A 1 = ±0.1 for large phases 0i ~ ±2. 
These asymmetries could be observable with 1.8 standard deviations after two years 
of running with the proposed luminosity. It can also be seen that in similar scenarios 
where Xi and \\ have larger Higgsino components the asymmetries could be even larger, 
and observable as well. The results here can be compared with the CP asymmetry in 
X1X2 production [32] in the same SUSY scenario. This analysis will be presented 
elsewhere [33]. In X1X2 production, although the cross section is much higher, the 
asymmetry is smaller and there are large backgrounds that further reduce it. The 
maximum statistical significance obtained for the asymmetry is 1.5<r when 0i = 37r/4. 
At any rate these two processes are complementary, because the asymmetries have 
a different dependence on 0i in each case. For instance, the CP asymmetry in X1X2 
production almost vanishes for </>i = n/2 and 0i = 37r/2 [33], while in selectron cascade 
decays it is nearly maximal. 

To conclude, we stress that the existence of CP-violating phases in SUSY-breaking 
terms of the Lagrangian is still an open question. In this paper we have focused on 
the phases 0i and M of the parameters Mi and /z, respectively. Assuming a relatively 
light SUSY spectrum there are two alternative possibilities to explain the unobserved 
EDM of the electron: either Mi and [i are approximately real or on the contrary 
they have large phases and there exist cancellations between neutralino and chargino 
contributions to the electron EDM. In this situation, the study of observables with a 
different functional dependence on 0i and as the CP asymmetry investigated here, 
is of great help in order to confirm one of the two hypotheses. 
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